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ABSTRACT 

We present if -band spectra of newly born OB stars in the obscured Galactic giant 
H II region W51A and ~ 0.8" angular resolution images in the J, H and Ks-bands. 
Four objects have been spectroscopically classified as O-type stars. The mean spectro- 
scopic parallax of the four stars gives a distance of 2.0 ± 0.3 kpc (error in the mean), 
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significantly smaller than the radio recombination line kinematic value of 5.5 kpc or the 
values derived from maser propermotion observations (6-8 kpc). The number of Lyman 
continuum photons from the contribution of all massive stars (NLyc ~ 1.5 xlO 50 s -1 ) 
is in good agreement with that inferred from radio recombination lines (NLyc = 1.3 
xlO 50 s _1 ) after accounting for the smaller distance derived here. 

We present analysis of archival high angular resolution images (NAOS CONICA at 
VLT and T-ReCS at Gemini) of the compact region W51 IRS2. The -f^-band images 
resolve the infrared source IRS 2 indicating that it is a very young compact HII region. 
Sources IRS2E was resolved into compact cluster (within 660 AU of projected distance) 
of 3 objects, but one of them is just bright extended emission. W51dl and W51d2 were 
identified with compact clusters of 3 objects (maybe 4 in the case of W51dl) each one. 
Although IRS 2E is the brightest source in the i^-band and at 12.6 fim, it is not clearly 
associated with a radio continuum source. Our spectrum of IRS 2E shows, similar to 
previous work, strong emission in Br7 and Hel, as well as three forbidden emission lines 
of Felll and emission lines of molecular hydrogen (H2) marking it as a massive young 
stellar object. 

Subject headings: HII regions — infrared: stars — stars: early-type — stars: funda- 
mental parameters — stars: formation 



1. Introduction 



The study of Giant HII regions (GHII - here taken to emit at least 10 LyC photons s , or 
~ 10 x Orion) in the near-infrared can address important astrophysical issues. The exploration of 
the stellar content o f obscured G a lactic GHII regio n s has been studied recently b y several groups: 



Han son et al 
fcoocl ) and lBik et al 



(2005. 



JOOjI). These observations revealed massive star clusters at the centers 



dl997h .Blum et al (jl999U200d . l200lMFigueredo et all (|2002l . l2005h . lOkumura et al 



of the HII regions which had been previously discovered and studied only at much longer radio 
wavelengths. 

The first step needed to trace the spiral structure of our Galaxy based on the spatial distribution 
of giant HII (GHII) regions is to derive spectroscopic parallaxes of OB stars in each embedded 
cluster. This is a well tested method in the optical window. We work in the near-infrared to 
overcome the obscuration at the optical wavelength s (A^/Ay 0.1) an d use the recently developed 



spectral classification scheme in the near infrared ([Hanson et al.lll996l ) to identify and classify the 
specific types of massive stars present. The uncertainty in the dist ance estimate is d ominated 
by the scatter in the intrinsic My of O-type stars (± 0.67 mag - IVacca et al.lll996l ). leading 



to an uncertainty of 30% for a single O-type star. The uncertainty is diminished as spectra are 
accumulated, and can reach <10% if enough stars are observed. If B-type stars are included the 
individual accuracy is better as their Mk are more constrained given their longer main sequence 
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lifetimes. 



Our procedure is independent of the traditional kinematic distance method based on radial 
velocity measurements of radio recombination lines and adoption of a Galactic rotation model. 
Almost all of the currently determined distances of Galactic GHII are based on this "kinematic" 
procedure which contains an ambiguity for directions inward of the solar circle and uncertainties 
introduced by random velocity c omponents in the gas. In our previous work we were able to break 
the distance ambiguity of W31 (|Blum et al.l 120011). placing it o n the near side of the GC at 3.4 
± 0.5 kpc from the Sun. In the case of W43 (jBlum et alJll999h and W42 ([Blum et al-lboool ) we 
also derived smaller di stances than those inferred from the pure rotational model; in G333. 1-0.4 
(jFigueredo et al.l 120051 ) , we also derive a smaller distance than given by the kinematic method 
(2.6 ± 0.4 kpc). Similar resul ts have been found by other investigators using the same technique 
(jHanson et al.lll997l : lBildl2004l ) for a number of other GHII. A recent H? Q maser study of the star 
forming HII region W3 finds a similar result based on maser parallax (jXu et al.ll2006l ): i.e., W3 
is closer than the kinematic distance and consistent with OB star absolute magnitudes in nearby 
associations lending credence to our results. 

W51 is a co mplex of giant HI I regions that was first detected by its free-free emission in the 
radio continuum (|Westerhoutlll958l ). The radio observations show that W51 is, in fa ct, composed o f 
a group of individual sources spread in an area of one degree in the Galactic plane (jBiegingl Il975l ) . 
The most luminous region in the W51 complex is G49.5-0.4 (R.A.= 19h23m42.02 and Dec. = 
+14d30m33.56s J2000), which is divided into eight smaller radio sources: W51A to W51H, ordered 
by R.A. W51A, one of the most luminous regions in the W51 complex , is lo cated at a kinematic 
distance 5.5 kpc (near distance), adopting the value given by iRusseill (|2003l ) who used the 1985 
IAU recommended values for the Sun to Galactic Center (GC) distance and rotation speed of the 
Milky Way, which are Ro = 8.5 kpc and 8q = 220 km s" 1 . Although it is one of the most luminous 
star-forming regions in the Galaxy, this source is completely obscured in the optical wavelengths 
due to the large extinction at optical wavelengths along the line of sight. 



Goldader &: Wynn- Williams! (|1994j) have imaged an area of 90" x 90" in the direction of W51A, 
using an angular resolution of ~ 0.4" in the near infrared. In addition to the radio sources that had 
been previously discovered, IRS1 and IRS2, they identified an third bright source - IRS3. They 
also argued that IRS2 is in fact a small cluster with at least twelve point sources in .fT-band. They 
have compared the images from a narrow band filter, centered in Br"f, with .ff-band images, and 
they have concluded that the extended inter-cluster emission is thermal from ionized hydrogen and 
neutral helium. The surface brightness of the extended emission in the radio and near-infrared 
led to a extinction of Ak s = 2.6 ± 0.3 mag for W51A, with areas still more obscured. We used 
this work for the identification of candidates for spectroscopy in order to determine the distance of 
W51A by the spectroscopic parallax method. 



Okumura et al.l (|2000l ) studied the properties of the stellar content of W51A in order to deter- 



mine the stellar masses, ages, spatial distribution and IMF. They used J, H and K imaging on an 
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area of 15' x 15' covering the G49.5-0.4 complex. These authors derived the interstellar extinction 
in the W51A direction by using the extinction in the visual band determined in a control field 
20' from the source IRS2 along the Galactic plane. The histogr am of the numbe r of st ars in each 
extinction bin shows two peaks at Ay = 7 and Ay = 19 mag. lOkumura et al.l (|2000l ) explained 
this effect as the presence of the Sagitta rius and Perseus spiral arms . By consideri ng the distances 



of five kpc and 13 kpc for the two arms (jGeorgelin &: Georgelinlll976l ). respectively, lOkumura et al 



( 20001 ) derived an aver a ge int erstellar differential extinction of 1.2 and 1.5 mag kpc -1 in the visual 
band. lOkumura et al.1 (120001 ) also found that the IMF slope in W51A agrees with the IMF slope 
in the control field in the range of masses between 10 and 30 M . On the other hand, they found 
an excess of stars more massive than 30 Mq. 

In the present work, we are interested primarily in constraining the distance to W51A using 
spectroscopic parallaxes of n ewly born OB stars. In fact, with 8.71 x 10 50 Lyman continuum 
photons per second (=NLyc 



Conti fc Crowther 



stars or ~ 200 07V-type stars ([Martins et al 



2004! ). this would be equivalent to ps 20 03V-type 
We present an investigation of the stellar 



content of W51A through JHK$ imaging and .fT-band spectroscopy (described in §2). In §3 we 
consider the photometry results, and in §4 we analyze the spectra. We determine the distance to 
W51A in §5 and discuss the results in §6. 



2. Observations and Data Reduction 

if -band (A ps 2.2 /jm, AA ps 0.33 //m) spectra of five objects in W51A were obtained on the 
nights of 2004 May 08 and 09 and 2004 June 01 using NIR0 on Gemini-North in queue-scheduled 
mode (GN-2004A-Q13). NIRI in f/6 mode delivers a plate scale of 0.12" pix" 1 . The f/6 camera 
and a four-pixel long slit (0.48" x 110") were used. The K grism produces if -band spectra with a 
resolution of R ~ 780 for NIRI and the linear dispersion is X/pix ps 7.08 A/pixel. 

J (A ps 1.25 fim, AA ps 0.2 fim), H (A ps 1.64 fim, AA ps 0.3 //m) and K s (A ps 2.15 /jm, AA 
ps 0.3 /xm) images of W51A were obtained on the night of 2005 May 07 usi ng the facility infrared 
camera ISp{] at the CTIO Blanco 4-m telescope (see brobst et all ((2003)). ISPI in F/6 mode 



delivers a plate scale of 0.3" pix 1 and a field of view of 10.25' x 10.25'. 

All basic data reduction was accomplished using IRAF@. Each image was flat-fielded using 
dome flats and then sky subtracted using a median-combined image of independent sky frames 
obtained 10' North and 10' East from W51A for direct imaging and three dithered frames l'-2' 



Gemini's Near InfraR ed Imager and spectrograph was built by the University of Hawaii's Institute for Astronomy. 
See iHodapp et al I |2003l ) for more details about the instrument. 

2 ISPI (Infrared Side Port Imager - "eye-spy"). More details about the instrument are available at the CTIO 
webpage. 

3 IRAF is distributed by the National Optical Astronomy Observatories. 
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far from the target for spectroscopy. One dimensional spectra were obtained by extracting and 
summing the flux in ± 4 pixel apertures. The extractions include local background subtraction 
from apertures ~ 1" on either side of the object. 



2.1. Photometry 

The images of W51A were obtained in non-photometric (cirrus) conditions. Total exposure 
times were 15, 3.4 and 2.4 minutes at J, H and Ks, respectively. The individual J, H and K$ 
frames were shifted and combined. These combined frames have point sources with FWHM of ~ 
0.92", 0.91" and 0.78" at J, H and K s , respectively. 

Figure Q] shows a false color image, made by combining the three near infrared images and 
adopting the colors blue, green and red, for J, H, and Ks, respectively. This Figure shows two 
regions with strong Br7 emiss ion probably associat ed with embedded star clusters seen in the 
center (W51A - regi on 3 from lOkumura et al,l |2000| ) and in the NE of the image (region 2 from 



Okumura et al.ll2000l ). The bluest stars are likely foreground objects, and the reddest stars are 
probably ifs-band excess emission objects, indicating the presence of hot dust associated with 
recently formed stars in the cluster. Background objects seen through a high column of interstellar 
dust would also appear red (but also fainter). The rec tangle in Figure [Q indicates the region, 
including W51A, where DoPHOT (jSchechter et al.lll993l ) photometry was performed on the com- 
bined images. DoPHOT detected a total of 2419 objects inside this rectangle in the Kg-band 
image. Of those 2419 objects, 1338 were also detected in ff-band and 855 in all three filters. All 
objects detected in J or ff-bands that were not picked up in K were excluded from our catalog 
because their errors were bigger than the cutoff limit (see below) . Figure [2] shows a finding chart 
using the -Kg-band image corresponding to the rectangular area of Figure HJ The numbers in the 
figure correspond to the brighter sources detected by DoPHOT. 

Standard stars were not observed but flux calibration was accomplished using stars in the same 
field, available in the 2MASS9 All-Sky Point Source Catalog (PSC). Twenty uncrowded stars on 
the W51A images were used for this purpose. The DoPHOT magnitudes were compared to the 
2MASS apparent magnitudes in order to define a zero point to the photometric calibration. The 
rms scatter in the zero point (difference between DoPHOT and 2MASS) in J, H, and Ks was 
0.032, 0.051, and 0.082, respectively. The error in the mean for the twenty stars is then 0.007, 
0.011, and 0.019 in J, H, and Ks, respectively. Uncertainties for the J, H and Ks magnitudes in 
W51A include the formal DoPHOT error added in quadrature to the error in the mean of the zero 
point offset. The sum in quadrature of both uncertainties results in typical bright star uncertainties 
of ± 0.019, ±0.028 and ±0.038 mag in J, H and Ks, respectively. The derived magnitudes and 



4 2MASS: Two Micron All Sky Survey 
http: / /www. ipac.caltech.edu / 2mass /releases/allsky / 



Skrutskie et al 



( 200rJ ) 



available 



at 
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colors of the stars labeled in Figure [2] are given in the Table [TJ 

The completeness of the DoPHOT detections was explored through artificial star experiments 
in the rectangular region of Figure 1 (see the discussion in Figueredo et al. 2005 for a description 
of the technique). The completeness is better than 90% for J and H < 16.5 and for K$ < 16. In 
other words, the detection of faint sources is more complete in J and ff-bands than in K$- This 
is likely because we are probably accessing mainly the (dense) foreground population and most of 
the fainter cluster members are not detected at this angular resolution (i.e., the true completeness 
limits in the clustered regions of the image are likely much brighter). The strong nebulosity will 
also hamper the detection of faint objects in the cluster region. Our primary goal in this paper is 
to derive the distance to W51A, and hence we only rely on the much brighter massive members of 
the cluster where crowding is not a significant factor. Indeed, the spectroscopic targets described 
below are in less crowded regions of W51A. 

The mean of the differences between the input and output H and K$ magnitudes of the 
artificial stars experiment, considering only the range of magnitudes where the completeness is 
better than 90 %, are ~ 0.07 and . 11, re spectively. The results of this experiment are similar to 



the ones showed in iFigueredo et al.l (|2005l ) . The magnitudes errors determined by this method are 
much more realistic than the formal errors defined by DoPHOT. We adopted an arbitrary cutoff 
of 0.2 mag (stars with larger errors were excluded from further analysis). 



2.2. Spectroscopy 

In order to avoid spectral degradation by bad pixels, the total exposure time needed to obtain 
S/N ~ 100 was split and individual spectra were obtained at seven different positions along the 
slit. After each exposure, the telescope was offset by three arcsec along the slit. In two cases, two 
or three stars could fit in the same slit, therefore the exposure time was adjusted to the fainter star. 
We needed blank sky spectra in order to subtract the sky emission of each target spectrum that is 
contaminated by nebular emission. The sky background were obtained by median combination of 
three dithered exposures in a blank field one or two arcminutes from each target (total = 10 sky 
positions). Following this procedure, the K-bsnd spectra of five of the brightest stars in the W51A 
clust er were obtained: RS6 (in this pape r: #50 - see Figure |2j), RS7 (#44) and IR S2East (IRS2E) 



(200 



from Goldader fc Wvnn-Williamsl (jl994h and #2519 (#61) and #2670 (#57) from lOkumura et al 



One dimensional spectra were obtained by extracting and summing the flux in ± eight pixel 
apertures (0.96") centered on the target in each dithered position. The extractions include local 
background subtraction from apertures ~ 1" on either side of the object. Wavelength calibration 
was accomplished by using the spectrum of an Argon lamp obtained each night. The error of this 
calibration is < 0.8 A (~ 10% of a pixel). In order to remove telluric features in the spectra of W51A 
stars, the A0V stars HIP92177 and HIP98640 were observed using an identical procedure, just after 
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or before those obtained for the target. At the time of the observation, both AOV stars have an air 
mass very close to the mean air mass attained by the targets in W51 during the observation period 
(Delta (A) ~ 0.02 airmass). The Br7 photospheric feature was removed from the average AOV-type 
star spectrum by fitting and subtracting a Lorentz profile. Since the Br7 region is relatively free 
from strong telluric features, this procedure is enough to obtain the template for telluric lines. The 
wavelength calibrated spectra were then divided by the continuum of the AOV-type star to remove 
telluric absorption. 

The JT-band classification scheme for OB stars is based on the (relatively weak) lines of CIV, 
Hel, NIII and Hell. The spectra of young stars in HII regions are often contaminated by the 
2.058 //m Hel and B r7 nebular lines but these features are not used for classifying O-type stars 
(jHanson et al.lll996l ). 



After the observations were obtained, it was found that lMader &; Angiond (|1996l ) have identi- 
fied the "AOV" telluric star HIP92177 (KOAql) as a binary system, and the secondary star is cooler 
than a GO-type star. Fortunately, the spectra of G-type stars don't show any strong features near 
the most important lines used to identify the O-type stars that are located between 2.05 fim and 
2.2 jjLva.. In fact, the Calcium, Iron and Sodium lines, t ogether with the CO absorption bands, 
common in cool stars (ILancon Sz Rocca-Volmerangelll992l ). are located at wavelengths longer than 
2.2 fj,m. In this work we thus consider the star HIP92177 as a telluric star, ignoring all lines longer 
than 2.2 fim wavelength. 



3. Results of the Photometry 



The ISPI J, H and -fTg- band images reveal W51A as an embedded star cluster readily seen in 
the center of Figure O Due to the strong nebulosity it is hard to detect the faintest stars in this 
region. On the other hand, a comparison between the density of stars in the center of this figure 
with a sky image close to W51A, allow us to verify that the stellar density is higher in the center 
than the other areas of the image. 



Okumura et al.l (|2000l ) found W51A is younger than 0.8 million years. This is in agreement 
with the photometric characteristics of this region, such as the high local reddening indicating that 
the massive stars have not completely shed their birth cocoons, and it is likely they are on the zero 
age main sequence (ZAMS). The cluster reddening ca n be approximate d as Ak ? ~ 1-7 x (H — 



since the intr insic colors of OB stars are nearly zero (|Koornneea 1 19831 ) ; see ICardelli et al.l (jl989l ) 
Mathisl jl99ch . 



The H — Ks versus K$ color magnitude diagram (CMD) is displayed in Figure El The labels 
in all plots refer to the same stars as in Figure [2j Both vertical lines indicate the position of the 
theoretical ZAMS shifted to 2.0 kpc (§6) and with interstellar reddening Ak s = 0.65 mag (dashed 
line) and an additional "intra cluster" reddening component of Ak s = 3.27 mag {Ax s totai = 
3.92 mag - solid line). The strong nebulosity in the cluster region, together with the crowding, can 
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explain the gap at H — K$ ~ 2 that appear in Figure [3] between K = 13 and K 
This gap is probably generated by the incompleteness as described in §2.1. 



15 magnitudes. 



The J — H versus H — Kg color-color plot is displayed in Figure HI In that diagr am the solid 



lines , from top to bottom, indica te interstellar reddening for main sequence M-type (jFrogel et al 



19781 ). O-type (jKoornneedll983l ) and T Tauri (IMever et all 1 19971 ) stars (dashed line). The solid 
vertical line between the main sequence M-type and O-type lines indicates the position of ZAMS. 
Asterisks indicate Ak s = 0, 1, 2 and 3 extinction values (in magnitude). Dots are objects detected 
in all three filters. The error bars in both figures refer to the final errors in the magnitudes and 
colors. The small number of objects between both solid lines in Figure H] means that most of 
stars have not been detected in the J-band images due the high extinction at this wavelength 
(Aj ~ 10.5 mag). Objects found to the right of the O-type star reddening line in the CMD of 
Figure HI shown as a solid diagonal line, have colors deviating from pure interstellar reddening. This 
is frequently seen in young star clusters and is explained by hot dust in the immediate circumstellar 
environment. 

For cases where the object has an apparent J — H less than an unreddened O star, the 
explanation is not clear. Most of cases, but not all, those objects are faint stars. Such objects may 
actually be blends of two nearby stars of different colors which are inadvertently matched when 
merging the photometry from different filters. This can happen when a red cluster star is projected 
near to a bluer foreground star. 



3.1. Reddening and Excess Emission 



Two concentrations of objects are seen in the CMD around the two vertical lines. The first 
one appears around H — K$ ~ 0.4, which corresponds to an ext i nction of Ak s = 0.65 mag 
(Ay ~ 6.5 mag) using the interstellar reddening curve of iMathisI (|l990l ). This concentration 
represents foreground stars (see also Figure H]). The foregroun d extinction can be estimated a s 
Ak s ~ 0.3 mag/kpc when using the distance of 2.0 kpc to W51A. iJonch-Sorensen &: Knudd (|1994l ) 
found smaller values, Ay « 1.8 — 2.5 mag /kpc or Ak ~ . 18 — 0.25 mag/kpc, for a field close 



to the galactic plane. On the other hand, lOkumura et al.1 (120001 ) defined as foreground sources 
in the direction of W51A cluster, all objects with Ak < 0.9 mag/kpc. The second and smaller 
concentration of objects appears around H — K$ ~ 2.0 or Ak s = 3.27 mag, indicating the average 
color of cluster members. A number of stars display much redder colors, especially brighter .fT-band 
sources. These objects are located at H — K$ > 2.0 and Ak s > 3.2. 

Determining the local reddening of a young cluster such as W51A is not simple, as the con- 
centration of stars does not appear as a clear main sequence, and the scatter about the mean value 
of Ak s is high. In fact, the differential "intra-cluster" reddening and the inclusion of very young 
sources produces the variation of the points on the CMD for H — K$ > 0.4. The local reddening can 
be determined through the colors of the O-type stars we have classified (see §4.1). The mean color 
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of the four stars spectroscopically observed are (H — Kg) = 2.25, which gives Ak s = 3.92 mag. 
The resulting local reddening in W51A direction, after subtracting the interstellar component, is 
then Ak s ~ 3.27 mag (Ay ~ 33 mag). 

In order to place the ZAMS in the CMD, as Figure [3] shows, the reddening and the distance 
of W51A must be derived. The procedure used to transform bolometric magnitudes (M& z) and 
effective temperat ures (T e ff) into colors {H — K) and absolute K magnitudes has been given by 



Blum et al 



Blum et al 



2000J). No transformation has been derived between ISPI and the system used by 
(|2000i ) (CIT/CTIO). The difference between both systems are negligible in comparison 
with the final uncertainties in the photometry. Figure [3] shows the ZAMS as a vertical dashed 
line, shifted to D = 2.0 kpc and reddened by Ak s = 0.65 mag due to the interstellar component. 
When adding the average local reddening (Ak s =3.27 mag), the ZAMS line is displaced to the 
right and downward, as showed by the solid line. A 5.5 kpc distance for W51A would move this 
ZAMS downward in the CMD by two magnitudes and the observed magnitudes for the spectroscopic 
targets would be inconsistent with a ZAMS luminosity class. 

Besides the strong local reddening in the direction of W51A, the color-color diagram and 
CMD show several objects whose colors deviate from the pure interstellar reddening. These objects 
are located to the right of the O star reddening line in the CMD of Figure El We can estimate 
a lower limit to the excess emission in the i^-band by supposing that the excesses at J and H 
are negligible, and that the intrinsic colors o f the embedded s t ars co rrespond to normal OB stars. 
Indeed, we can follow the same procedure as lFigueredo et al.1 (J2002J) by assuming that our sample 
of stars is composed of young objects (not contaminated by f oreground or bac kgro und stars) which 
have i ntrinsic colors in the range (H — K) Q = 0.0 ± 0.06 mag (|Koornneeflll983l ). In lFigueredo et al 



(|2005l ) we adopted for all objects the intrinsic colors of a B2 V star and from the difference between 
the observed J — H and the adopted B2 V intrinsic J — H color, we obtained the J-band extinction 



by using the adopted extinction law (jCardelli et al.l Il989l ; iMathisI Il990l ) . The results for the same 
analysis in W51A are displayed in Figure [5j 

In Figure [5] the solid line indicates K exc = 0. Connected solid diamonds refer to the average 
value of the K excess in 1 magnitude bins. Dashed lines indicate 1, 2 and 3 a from the average. 
Bright objects with very large excess (very positive values) in Figure [5] cannot be explained by 
errors in the dereddening procedure. The excess should be due to accretion disks and/or hot dusty 
envelopes around the less massive objects in the cluster. A few objects are well above the 3 a 
scatter for otherwise normal stars in Figure [5j 

As with in the case of G333. 1-0.4, most of the stars in W51A (FigureEJ) have a modest n e gative 



excess, about 0.2 magnitude. This negative excess has been explained by iFigueredo et al. 



as a consequence of our assumption that all stars have the intrinsic color of a B2 V type star; in 
reality, most of the stars will have slightly positive intrinsic color. Our goal is to identify stars with 
a significant excess which would cause them to lie to the right of reddening line in Figure [U so 
this modest negative excess for "normal" stars, or stars with a small excess is not important for 
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our purposes. In the following sections we have only corrected the if-excess for stars with positive 
excess as determined here, for all others we impose zero excess emission. 

It is important to note that the excess emission as determined in this section is only a lower 
limit, since we assumed the excess was primarily in the if- band (but according to F igure El there are 
not many stars with large excess for the higher masses). iHillenbrand et al.l (|1992l ) have computed 
disk reprocessing models which show the excess in J and H can also be large for disks which 
reprocess the central star radiation. The ifs-band excess e mission determined for IR S 2E, a bout 
two magnitudes, is in agreement with the values found by IHillenbrand fc Carpenter! (|2000l ) for 
young stars in the Orion cluster. Objects on the top of Figure [5] show excess emission still higher 
and it is an indication that probably such peculiar objects a re, in fact, massive YSOs or ultra- 
compact HII (UHII) regions, as should be the case of IRS2E. iGoldader &, Wynn- Williams! ()1994l ) 
have found, besides IRS2, a pair of luminous sources in W51A - IRS1 and IRS3. The coordinates of 
IRS1 match those of #403 in figure El and IRS3 with #96 and #152, probably IRS3E and IRS3W 
respectively. Only #152 has been detected in H and if-band. Lower mass protostars also can 
show excess emission in .fT-band. 

As we can see in Figures [31 [H and [5l some objects in W51A are very bright, display very red 
colors and at the same time have excess emission well above the 3 a scatter for the average of 
normal stars. This is the case for objects such as IRS2W. 

Given suc h evidence for circ u mstel lar disks, we attempt to estimate by iteration the excess 
emission using IHillenbrand et al.l (119921) m odels for reprocessing disks, such as was done in the 



case of G333. 1-0.4 (jFigueredo et al.l 120051 ). to get a sense for how bigger the excess might be 
compared to that derived by assuming all the excess is in the Ks-band. This was done for the 
star we have detected in all three bands, IRS2W. The results corresponds to an excess emission of 
K e xc = 3.58 mag from a face-on reprocessing disk with a central source of 15 M & corres ponding 



to the spectral type BO. 5. The resultant reddening for this star was determined by using iMathie 
(|l99oh approximation as A Ks « 1.7 x ({H - K s ) + 0.04 - 0.5). The resultant extinction for IRS2W 
after correcting for the excess as determined above, is Ak s = 4.25 mag. 



4. Results: Analysis of Spectra 



The Gemini North/NIRI spectra of five objects in W51 sources are shown in Figures [6] and [71 
All sources have been divided by a low-order fit to the continuum after correction for telluric ab- 
sorption. The signal-to-noise ratio is S/N « 100. These spectra have been background subtracted 
with nearby (~ l".0) apertures, though non-uniform extended emission could affect the resulting 
He I and Br7 absorption apparently present in the stars themselves. 

In order to identify the spectral t ype of the sta r s we c ompa red all spe c tra w ith the iC-band 
spectroscopic standards presented by lHanson et al.l (|1996l ) and iBik et al.l (|2005l ). The features 
of greatest importance for classification are (vacuum wavelengths) the CIV triplet at 2.0705 fim, 
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2.0769 /mi and 2.0842 /im (emission), the NIII complex at 2.116 /im (emission), and Hell at 2.1891 
/im (absorption). Stars cooler than spectral type 07 show Hel (absorption) to the blue side of the 
NIII complex, at 2.1126 /im (which can be confused with nebular over-subtraction). 



The NIII complex has been identified by Hanson et al.l (119961 ) in early-0 stars as being due 
to NIII(8— 7). The fact that this transition is a multiplet system explains the broad profile of the 
feature, centered at 2.116 /im. It shows constant width through moderate variations in temperature 
and large variations in luminosity. The CIV triplet (2.071 /im, 2.077 /im and 2.084 /im) is typically 
weak and seen only in very high signal-t o -noise spectra ([Hanson et al.lll9 96). The present classi- 
fication system laid out bv lHanson et al,l (|1996l ) does not have strong luminosity-class indicators. 
Still, without nebular contamination (which is not the case presently), the Hel (2.0581 /im) and 
Br7 (2.1661 /im) features can be used to app roximately distingu ish between dwarfs plus giants on 
the one hand, and supergiants on the other ([Hanson et alJl2005l ). Generally strong absorption in 
Br7 is expected for dwarfs and giant stars and weak absorption or emission for supergiants. 



4.1. O type Stars 

Figure [6] shows four sp e ctra t hat may be compared with the -fT-band spectroscopic standards 



presented by lHanson et al.l (|l996l ). The signal-to-noise ratio is S/R ~ 100 — 120, high enough to 
separate and identify the CIV triplet and NIII complex lines. The telluric absorption lines have 
been removed from the #44, #50 and #61 original spectra by using the HIP92177 spectrum. The 
spectrum of star #57 has been divided by HIP98640. The measured equivalent widths for the main 
lines for classification of the stars in Figure [6] are shown in the Table [2j The classification resulted 
in the spectral types indicated in the Figure. 

The CIV triplet at 2.069-2.083 /im (emission) and the NIII complex at 2.116 /im (emission) 
identify the four sources in Figure [6] as early to mid O-type stars. The relatively strong NIII and 
Hell lines that appear in the spectrum of #57 indicate that this star is hotter tha n 06. On the 



other hand, the absence of CIV at 2.071 /im places #57 in the k03-04 subclass ([Hanson et al 



19961). The CIV line appears very weak and blended with the third line of the triplet. The CIV line 
in the spectrum of the source #44 indicates that this star is cooler than a typical 04 and hotter 
than an 07-type star. The NIII line is still stronger than in the previous star, placing #44 in the 
k05 subclass. In the spectrum of #50, the NIII line appears weaker than in #44 and for this reason 
this star is probably cooler than the previous one. The CIV triplet is still present, but weaker than 
in #44, placing #50 between the spectral type 06 and 07. The combination of features for source 
#61 suggest it is the coolest one in Figure [H and we place it in the k07-08 subclass. As in our 
earlier work, we adopt a ZAMS classification due to the presence of massive YSOs in the cluster. 
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W51 IRS 2 



5.1. High Resolution Images 

W51 IRS 2 (hereafter IRS 2) consists of a number of a stellar objects surrounded by a peanut- 
shaped extended emissio n (as seen in low resolution images at 2.2 //m) of ~ 15" x 15". Two bright 
sources were reported bv lGoldader &; Wvnn- Williams! ()1994l ) from low resolution near-infr ared im- 
ages: IRS 2E and IRS2 West (hereafter IRS 2W). IRS2 is also a radio sourc e, named W51d (IMartin 
19721 ). however the radio peak does not coincide with the infrared peak (jGenzel et al.l ll982) IRS 
2E. 

In this section, we present high resolutio n images of IRS2 t a ken with the adaptiv e optics near 
infrared camera NAOS CONICA (NACO) (jLenzen et al.l I2OO3I : Eousset et alJ bood ) at the ESO 
VLT UT4. The raw data were retrieved from the public database archive (ID 71C-0344(A)). We 
processed the data combining the dithered images and subtracting sky images taken with the same 
exposure time, but displaced from the target. The images were flat-fielded and a bad pixel map 
was applied to th e final combined im ages. Photometry was performed by means of the IDL code 
XSTARFINDER (jDiolaiti et al.ll2000l ). a psf fitting algorithm suitable for adaptive optics assisted 



obser vations. The extracted fluxes were calibrated using the standard star P9150 (IPersson et al 



19981 ) observed as part of the nightime calibration set. The observing log of the images is presented 
in the Table O The zoomed box in Figure [T] is a false color image of IRS2, coded as red (K s ), green 
(H) and blue (J) and constitutes the highest r esolution image o btained from this region until now 
(the //—band image was presented recently bv lLacv et al.ll2007l ). 



The H and K images of IRS 2 are shown in the Figure El In both images, the circles represent 
the position of objects detected bv lGoldader &, Wynn- Williams! ()1994l ) through a PSF fitting algo- 
rithm. The position of IRS2E is marked by a box and, as can be seen, three sources lie within it. 
Goldader & Wynn- Williams concluded that IRS2E is a small cluster of stars. The high resolution 
images suggest a somewhat different picture. Three sources are clearly identified both at H and 
Ks- Two of the sources (a and c) appear stellar, while the third (6) is an extended clump of 
emission. Furthermore, source c appears relatively blue and may be a foreground object. Source a 
is very red and the brightest Ks source in the IRS2 region. We thus associate IRS2E with a single 
unresolved source (a). 

We present the IRS2 sources as open squares on the CMD in the Figure [3l We detected 36 
sources in all three bands, however due to the low quality fitting of the two brightest sources in the 
J-band, they were removed. The source IRS2E, on the other hand is not detected in our J-band 
image. 
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5.2. Radio Sources and MIR Emission 



IRS 2 hosts three UCHII regions (jGaume et al.lll99 3). The high resolution images show that a 
cluster of three objects can be associated with IRS2E within 1" (which w ould be 0,01 pc at a distanc e 
of 2.0 kpc). The main radio source W51d is a cometary UCHII region (IWood ChurchwellHlQ SQ). 
and th e other two sources are W51dl and d2 as showed in Figure El iGoldader fc Wynn- Williams 
(| 19941 ) noted that instead o f extended emis s ion, n o NIR counterpart was found at the radio peak 
position, on the other hand. lOkamoto et al.l (|200ll ) assumed that IRS 2W is the NIR counterpart of 
W51d, only by its proximity and by the visual inspection of the morphology of the extended MIR 
emission. Figure [9] is a composite image showing both radio and MIR emission overplotted on the 
Kg-hand image of IRS2. The heavy black contour lines represent MIR emission at arbitrary units, 
from high resolution images taken recently at the Gemini South Observatory at 12.6 /im by us. 
A MIR study of W51A will be presented in a forthcoming paper (Barbosa et al. in preparation). 
The light white c ontour lines represent the r adio emission, also at arbitrary units, taken from the 
radio catalogue of lWood &: Churchwelll (|1989l ). Taking the .ffs-band image as reference, we shifted 
the radio image on 0.8 arcsec to match the position of W51d and its NIR counterpart. This value 
is within the uncertainties of 1 arcsec as expected by UT4 pointing. The MIR image was shifted 
to match the MIR infrared peak of IRS2E and its NIR counterpart. The shift between MIR and 
NIR images amounts to 2 arcsec. The MIR images seem to have an uncertainty worst than the 
quoted 1 arcsec, specially at 25 //m, so we can consider this shift value is within the uncertainties. 
The position of W51d2 was derived from its r elative position to W 51di comparing the coordinates 
presented by IWood fc Churchwelll (|1989l ) and iGaume et al.l (|1993l ) for the former source. 



The composite image sh ows IRS2W as the NI R counterpart of W51d, but also shows that 
IRS2W is not a MIR source. lOkamoto et al.1 (|200ll ). on the other hand identified IRS2W as the 
NIR counterpart of their MIR source OKYM3. The Western MIR emission peaks at the position of 
two stellar objects that are also coincident with the secondary peak of radio emission from W51d. 
This peak is spatially resolved in the MIR, however, and it is not clear if there is a buried source 
associated with it or if it is associated with an ionized cavity produced by IRS2W. We have recently 
obtained a Gemini NIFS spectrum of IRS2W which indicates it is indeed an O-type star of the 
hottest spectral type. This spectrum will be the subject of a forthcoming paper. 

IRS2E is the brightest source in the i^s-band and at 12.3 /im as well, and it is not a radio 
source. This lack of emission at longer wavelengths may be due to its youth, it may be the 
case that this object has a low radio luminosity and/or a small emitting volume giving a low 



emissi on measure. However, it is worth noting that the radio maps presented by IGaume et al 



(jl993l ) show some extended emission near, but not coincident, with the position of IRS2E. Thus we 
conclude IRS2E is a massive young stellar object (MYSO), but we can not identify it as an UCHII. 
W51di is seen as a r a dio, N IR, and MIR source and W51d2 was not detected in the radio maps of 
Wood &: Churchwelll (|l989l ) and may have weak MIR emission. 



The high resolution images show that a cluster of three objects can be associated with IRS2E 
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within 1", which represents ~0.01 pc at the new distance of 2.0 kpc. It is also interesting to note 
that at the position of the radio source W51di a small cluster of at least four sources is seen. A 
closer inspection of the images reveals two other sources, but they are too faint to be detected 
by a PSF fitting algorithm. We conclude that a small cluster of 4, maybe 6, sources (all of them 
within 0.3" or 660 AU) is the NIR counterpart of the radio source W51 di. A similar conclusion 
can be taken from the near infrared counterpart of source W51 d.2, since three sources are seen in 
its position, all of them separated by less than 0.3". 



5.3. IRS2E Spectroscopy 



The spectrum of the most luminous source in the K$ band in our sample, IRS 2E, is shown 
in Figure [71 The position of this source in Figure [3] shows that IRS 2E is more luminous in K$ 
than the O-type stars on the ZAMS which we classified in the previous section. OB stars can 
be very luminous in the NIR during their stage as a MYSO. This high luminosity is due to the 
excess emission in the K band (Figure [5]) . This -KT-band excess has to be accounted for if one wants 
to estimate the spectral type of the MYSO using its luminosity. The strongest emission lines are 
Br7 (A =2.166 /im) and Hel (2.058 and 2.113 ^m). Figure [7J shows a series of H2 and [Felll] 
lines. Also, we confir m the detection of the "unidentifi ed" infrared line at 2.288 /im reported by 
Okumura et al.l (|200ll ). This line was discussed earlier by lGeballe et al.1 (|199ll ) who did not detect it 
in IRS2 (apparently due to their large aperture of 5") but found it w as common in planeta ry nebulae 



with middle excitation temperatures. Based on their observations. iGeballe et al. (119911) we r e able 



to estimate the ionization potential of the species giving rise to the line which led iDinersteinl ([200 ll ) 
to iden tify it with [Se IV] • T he line is now known to be commo n in the spectra o f planetary nebulae 
(PNe) ([Sterling et al.l 120071 ) . A second line was identified by IDinersteinl (|200ll ) near 2.199 //m as 
being associated with [Kr III]. This line can be confused with an H2 transition in lower resolution 
spectra. Our spectrum of IRS2E shows a line at this positi on and it is closer to the [K r III] position 
than the H2 position. Independently of the present work, iBlum &: McGregor! (120071 ) identified the 
[Se IV] (and [Kr III]) line in two UCHII regions. Thus the "unidentified" line(s) appears more 
common than previously thought and may be produced in higher excitation regions very close to 
newly formed massive stars. 

The spe ctrum of IRS 2E is very similar to the spectrum of the UCHII region G25.2-l.74 
presented by iBik et al.l (|2005l ). Given the lack of a readily identifiable radio counterpart, this 
may simply reflect the nebular nature of IRS 2E and the environment within which it is found. 
Disentangling the internal and external excitation of IRS2E is beyond the scope of this paper. 
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The Distance to W51A 



In the previous section we classified the spectra of five objects in W51A: four O-type stars 
and a MYSO. In this section, the distance to W51A is determined by using the spectroscopic and 
photometric results found for the four O-type stars. Table 0] summarizes the main properties of 
the O-type stars. We compute distances assuming the O-type stars shown in Figure [6] are zero-age 
main-sequence (Dzams - column 5 of Tabled]) or in the dwa rf luminosity class (co lumn 6). For the 
dwarf case , the d istance is determined using the My given b vlMartins et al.l (120051 ) and V — K from 



Koornneefi (|l983l ). For the ZAMS case, we use the Ma- fromlBlum et al.l (120001) and adjus t the M k 



to account for the difference in My between the IVacca et al 



199fih and rMartins et all <j200fJ ) O 



star properties. iMartins et al.l (|2005l ) have revised the dwarf O star My, and we have assumed the 
ZAMS My will be similarly shifted. For the spectral types considered here, the maximum change 
is to make the My 0.25 mag fainter. 

The distance estimates are shown in Table[4l For the derived spectral types, we obtain distances 
of 2.0 ± 0.3 and 2.4 ± 0.4 kpc for the ZAMS and dwarf cases, respectively. The lack of luminosity 
class indicators do not allow us to unambigously determine the distance, but the ZAMS value is 
preferred given the presence of massive YSOs in the cluster. The uncertainty quoted is the error 
in the mean. 

An extra uncertainty on the derived distance comes from the inte rstellar ext inction law. In 



this p aper, we adopted the in t erstel lar reddening curve of iMathia (Il990i ). However. iMessineo et al 
(|2005l ) and iNishivama et al.1 (120061 ) predict Ak s /Eh~ks ~ 1-5 m contrast with the value 1.7 



suge ested bv lMathisI (|l99oh . Higher values of interstellar extinction will result in a underestimate 
of the distance. The distances quoted above would be a factor of 1.28 larger for the extinction law 
predictiong a smaller Ak- 



6.1. The Lyman Continuum Flux 



Conti &; Crowtherl (|2004l ) determined the number of Lyman continuum photons emitted per 
second by W51A as 8.71 x 10 50 with the clust er located at a distan ce of 5.5 kpc. This corresponds 
to w 20 03V- type stars or « 200 07V stars (jMartins et al.1 120051 ). The new spectrophotometric 
distance we have determined by using four O-type stars is D = 2.0 ± 0.3 kpc, considerably smaller 
than 5.5 kpc. The new distance brings the number of Lyman continuum photons down to NLyC = 
1.15 x 10 50 . 



The Lyman continuum flux comes primarily from the four brightest stars which are well above 
our photometric completeness limit. We have calculated that the number of Lyman continuum 
photons from the contribution of all massive stars we have identified in the cluster is NLyC = 
1.5 x 10 50 photons s _1 . 



The observed NLyc photons derived from radio techniques is only a lower limit to that emitted 
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by the stars, since some may be destroyed by dust grains or leaked through directions of low optical 
depth. On the other hand, due to the high extinction and crowding in the region some ionising 
stars may not have been detected in the J, H and ET-bands. Therefore, both values found for the 
total NLyC produced by the hot stars in W51A, by radio techniques and near-infrared, should 
be considered as lower limits. The NLyC photons from the radio measurements and that derived 
from the star counts are in good agreement. 



7. Discussion and Summary 



We have presented J, H and Kg images of the HII region W51A obtained on the 4-m Blanco 
Te lescope using the facili ty imager ISPI. The photometric results agree with the previous work, such 



as lOkumura et al.l (|2000l ). fT-band spectra of five of the brightest cluster members were obtained 
with Gemini North telescope using the instrument NIRI. Four of them have been classified as O— 
type stars and one as a massive young stellar object. We have derived the distance to W51A by 
using the spectroscopic parallaxes of four O-type stars. We also have presented archival images 
of the compact HII region IRS 2 taken with the adaptive optics camera NACO at the VLT. The 
high resolution images revealed a small cluster of four objects as the NIR counterpart of the UCHII 
region W51di and another cluster with three objects at the position of the UCHII region W51d2- 
The source IRS 2E was resolved into three sources and we identified the brightest source in the 
K band as IRS 2E. The K band spectrum of IRS 2E is very similar to that of an UCHII region, 
however, we do not explicitly classify it as so, because there is no distinct radio source with it. 

The value we have derived for the distance to W51A by spectroscopic pa rallaxes of O-type 
stars , 2.0 ± 0.4 kpc , is considerably smaller than the kinematic distance: 5.5 kpc (jConti fc Crowther 
20041 : iRusseill 120031 ) . The systematic uncertainties on the interstellar extinction law are not big 
enough to change this picture. As a consequence of this result, W51A is fainter than previously 
determined. Although the number of Lyman continuum photons is smaller with the new distance, 
~ 1.5 x 10 50 photons per second, W51A can still be classified as a giant HII region harboring a 
considerable number of massive stars. 

The newly derived distance i s much smaller than the d istance (7 ± 1.5 kpc) derived from H2O 
maser observations presented by lGenzel et al.l ()198ll . Il982l ) for the masers associated with W51A 
MAIN (ab out 1' south ea s t of I RS2E) and the similar re sult (8.5 ± 2. 5 kpc) for W51A NORTH 



derive d by ISchneps et al.l ()198ll ). The method used by iGenzel et al.l (119821') and ISchneps et al 



kpc| 
19821) 



()198ll ) assumes the maser motions are random and isotropic. Ilmai et al.l (|2002l ) present somewhat 
more detailed models of the maser kinematics in W51A NORTH (centered within 1" of IRS2E) 
and find that the motions are better explained by a radial out flow model which they use to fit the 
positions and outflow velocity. The distance derived from this technique is 6 ± 1.3 kpc, very close 
to the kinematic distance. On the other hand, they suggest the distance is completely covariant 
with the expan sion velocity and the direction of the flow and the uncertainty permits a distance 
value of 7 kpc. Ilmai et al.l (120021 ) further show that the masers for W51A MAIN and NORTH are 
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non-spherically symmetric and that there are significant residuals with respect to the assumption 
that all the maser spots emminate from a single point at the origin of the out flow (i.e. there are 
significant deviations for a simple radial out ward flow when looking at the derived maser space 
positio ns and velo c ity ve ctors). Using this assumption of random and isotropic motion for the water 



maser, 



Imai et al.l (j2002l ) estimated the distance of 4 kp c to W51A NORTH. This is in disagreement 



with the value found previously by lGenzel et al.1 ()198ll ). It is not clear what other models could fit 



the data or how serious the non-symmetry of the flow is in comparison to the model assumptions. 
In particular, it is not clear if a near distance model could fit the data. 

The situ ation may be even more complex in light of recent mid infrared observations (in [Ne 



II]) made by lLacv et al.l (|2007l ). These investigators discovered a large "fan like jet" or outflow 
which is ionized as it enters the HII region (apparenly from a source embedded in the molecular 
cloud) very close to W51A NORTH. They speculate that the ionized flow may have the same source 
as the W51A NORTH masers and this would require precession of the outflow since the ionized 
flow has a strong radial velocity while the masers are dominated by proper motions in the plane of 
the sky. 

The masers of W51A NORTH are probably associated with IRS 2, while the O stars are 
somewhat further away in projection from IRS 2. The masers and O stars may trace different 
physical locations along the spiral arm and appear close only in projection. Within the uncertainty 
of the maser techniques, this explanation is less favorable since the distance derived for W51A 
MAIN is similar to the larger distance to W51A NORTH and W51A MAIN is closer (in projection) 
to several of the O stars (see Figure [2]) . If the O stars are actually unresolved multiple (massive) 
stars they would appear brighter and thus apparently closer. On average, this would require that 
nine stars of similar luminoisty make up each object we have observed if we consider just the ratio 
of distances two kpc to six kpc. This seems an unlikely explanation. We are presently left with no 
satisfactory explanation for the distance discrepancy between the masers and radio observations on 
the one hand, and the near infrared spectrophotometric observations on the other. 

The authors wish to thank Ed Churchwell for making available the radio map of W51 IRS2 
in electronic form. We thank an anonymous referee for the careful reading of this paper and for 
the useful comments and suggestions which have resulted in a much improved version. AD and 
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Fig. 1. — False color image of W51A: J is blue, H is green and K$ is red. The red rectangle 
indicates the region where DoPHOT photometry was performed. The coordinates of the center of 
the image are RA (2000) = 19h23m42.02s and Dec. = +14d30m33.56s and the size of the image is 
10' x 10' (plate scale = 0.3" /pixel). North is up and East to the left. The zoomed box (top right) 
shows a VLT false color image of W51 IRS2 (see §5.1 and also Figures [8] and [9]) . 



-22- 




Fig. 2. — Finding chart using a i'Ts-band image of W51A, where DoPHOT photometry was per- 
formed (red rectangle in Figure [1]). Object labels refer to the brightest cluster members detected 
by DoPHOT. Squares refer to the masers W51A NORTH and MAIN (see §7). stars #44, #50, 
#57 and #61 whose spectra show photospheric lines typical of O-type stars, have been used to 
determine the spectroscopic parallax of W51A. The size of the image corresponds to an area of 
~ 5.1' x 3.5'. North is up and East to the left. 
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Fig. 3. — Ks vs H — Ks color-magnitude diagram (CMD) of W51A showing the position of the 
theoretical ZAMS shifted to 2.0 kpc and with interstellar extinction Ak s = 0.65 mag (dashed line). 
An additional "cluster" extinction component of Ak s = 3.27 mag (AK s totai = 3.92 mag) results in 
the ZAMS position indicated by the vertical solid line. Object number labels are the same as in 
Figure [2 Open squares indicate sources detected in IRS2 (see zoomed box in Figure [T] and §5.1 for 
more details). 



-24- 




Fig. 4. — J — H vs H — Ks color-color plot showing the reddening line of M-type stars {heavy 
solid line), O-type stars {solid line) and T Tauri stars {dashed line). Dots refer to stars detected 
by DoPHOT in all three bands J, H and Ks- The asterisks indicate the corresponding Ak s along 
the reddening vector. Object number labels are the same as in Figure [2] and Q] 
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Fig. 5. — Excess emission as a function of dereddened K$- band magnitude (K a ). Only objects 
detected in all three bands have been included. Circles indicate objects with measured J, H and 
Ks magnitudes. The solid line indicates K exc = 0. Connected solid diamonds refer to the average 
value of the K$ excess in one magnitude bins. A 2 sigma clipping was applied to calculate the 
average valuesDashed lines indicate one, two, and three a deviations from the average. IRS2W 
(Solid squares) have been detected in the VLT images (See §5.1 for more details). Very positive 
values represent circumstellar excess emission while small or slightly negative values do not; see 
text. 
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Fig. 6. — if -band spectra for the four brightest stars in the W51A: #57 (04V), #44 (05V), #50 
(06.5V) a nd #61 (07.5V). Th e spectral type was determined by comparison with the 0-type 



stars from lHanson et al.l (Il996l ). The spectral resolution of R ~ 780 gives a linear dispersion of 
~ 7.08 A/pixel. Spectra were summed in apertures 0".96 wide (eight pixels) by a slit width of 
0".48 x 110" and include background subtraction from apertures centered < 1".0 on either side 
of the object. Each spectrum has been normalized by a low-order fit to the continuum (after 
correction for telluric absorption). The spectra are often contaminated by the 2.058 fxva. Hel and 
Br7 nebular lines which can be over or under subtracted from the stellar source. The signal to 
noise ratio is S/N ~ 100 - 120. 
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Fig. 7. — if-band spectrum of IRS2E taken with the same setup as the O-type star spectra shown 
in Figure 
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Fig. 8. — W51 IRS2 in the if -band (left pannel) and Kg-band (right pannel). The circles labeled 
a, b and c refer to the position of IRS2E, an extended e mission (no stellar source) and a NIR 
source also detected by iGoldader fc Wvnn- Williams! ()1994l ) respectively. The star labeled "P" is 
also detected on the POSS red plate. The position of IRS2E is marked by a box (sized l"xl"). 
IRS2E is the brightest source in the field. The NIR counterparts of radio sources W51d, di and d.2 
are also indicated. 
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Fig. 9. — Compos ite image of IRS 2E . The i^g-band image is overplotted with radio emis- 
sion contours from Wood fc Churchwelll (|1989l ) in arbitrary units in light white contours from 



Wood &: Churchwelll (|l989l ) and by MIR emission contours (C. Barbosa, unpublished), also at ar- 
bitrary units in heavy black contours. Sources are indicated by black or white circles to give 
the best contrast. W51d 2 and the radio emission southern of W51dl were not detected by 



Wood fc Churchwell 



of iGaume et al.1 jl993 ) 



19891 ) . but fainter radio emission does extend toward IRS2E in the maps 



Table 1. Colors and Magnitudes* Catalog 





J 




H 




K S 




T — H 




IRS2E** 




> 21 


16.70 ±0.01 


10.80 ±0.01 


< 4.3 


5.90 ±0.01 


IRS2W** 


16.47 ±0.01 


14.66 ±0.01 


11.50 ±0.01 


1.81 ±0.01 


3.16 ±0.01 


#7 


17.09 


± 


0.02 


13.93 


± 


0.01 


11.92 


± 


0.02 


3.16 


2.01 


#17 


17.91 


± 


0.04 


14.36 


± 


0.02 


12.42 


± 


0.02 


3.55 


1.94 


#21 


17.13 


± 


0.03 


13.31 


± 


0.01 


11.33 


± 


0.02 


3.82 


1.98 


#24 


15.90 


± 


0.01 


13.16 


± 


0.01 


11.31 


± 


0.02 


2.74 


1.85 


#25 


17.79 


± 


0.04 


14.59 


± 


0.02 


12.67 


± 


0.02 


3.20 


1.92 


#32 


if? 17 
10. 1 1 


± 


n m 


i 9 no 


± 


n m 
U.U1 


1 O I'i 

iz. i 6 


± 


n no 
U.Uz 


2.25 


1.19 


#34 








16.31 


± 


0.13 


12.68 


± 


0.03 




3.62 


#44 


17.31 


± 


0.01 


13.67 


± 


0.01 


11.31 


± 


0.02 


3.65 


2.35 


#45 








15.15 


± 


0.02 


12.48 


± 


0.03 




2.67 


#50 


17.50 


± 


0.02 


14.11 


± 


0.02 


12.23 


± 


0.02 


3.39 


1.88 


#52 


13.67 


± 


0.01 


12.53 


± 


0.01 


12.09 


± 


0.02 


1.15 


0.44 


#57 


16.88 


± 


0.02 


12.91 


± 


0.01 


10.77 


± 


0.02 


3.97 


2.14 


#59 


13.24 


± 


0.01 


12.43 


± 


0.01 


12.16 


± 


0.02 


0.81 


0.27 


#61 








15.16 


± 


0.02 


12.51 


± 


0.02 




2.65 


#62 








14.20 


± 


0.02 


12.28 


± 


0.02 




1.93 


#64 


18.17 


± 


0.06 


14.36 


± 


0.02 


12.35 


± 


0.02 


3.80 


2.02 


#73 


16.27 


± 


0.02 


13.27 


± 


0.02 


11.57 


± 


0.02 


3.00 


1.70 


#98 








16.56 


± 


0.03 


13.00 


± 


0.02 




3.56 


#103 


16.96 


± 


0.01 


14.44 


± 


0.02 


12.83 


± 


0.02 


2.52 


1.61 


#145 








16.32 


± 


0.04 


13.16 


± 


0.04 




3.16 


#152 








16.86 


± 


0.06 


11.68 


± 


0.07 




5.18 


#754 








17.28 


± 


0.08 


12.94 


± 


0.08 




4.35 



*The sum in quadrature of the zero point uncertainty plus the PSF-fitting un- 
certainty define the final magnitude errors in J, H and K$, see §2. 



**VLT data 



Table 2. Equivalent width (A) of strategic lines used to classify the O-type stars 



Lines 


#57 


#44 


#50 


#61 


Hel (2.059/xm) 


1.75 ± 0.13 


1.42 


± 


0.22 


5.35 ± 0.37 


6.21 ± 0.28*ab 


CIV (2.071//m) 




0.37 


± 


0.05 


0.72 ± 0.14*b 


0.24 ± 0.12 


CIV (2.080Atm) 


0.32 ± 0.04*b 


0.96 


± 


0.08 


0.53 ± 0.14 


0.61 ± 0.17 


NIII (2.116/im) 


0.95 ± 0.05 


1.82 


± 


0.09 


0.64 ± 0.11*b 


1.70 ± 0.21 


Brg (2.166/um) 


0.68 ± 0.09 


0.41 


± 


0.09 


3.83 ± 0.22 


11.23 ± 0.20*ab 


Hell (2.189/im) 


0.21 ± 0.10*ab 








0.51 ± 0.16*ab 


0.63 ± 0.12*ab 



*b - blended lines 
*ab - absortion lines 



Table 3. NACO observations of W51 IRS2: instrumental set-up and observation logs. Camera 
S27 provides images with resolution of 27.15 mas/pixel and FOV of 28 x 28 arcsec employing an 

Aladdin InSb 1024 x 1024 pixel array. 



Filter 


Central Wavelength 


Bandwidth 


DIT a 


NDIT b 


PSF c 


Date of 


Seeing d 




(H 


(H 


(s) 


(s) 


(arcsec) 


observation 


(arcsec) 


J 


1.27 


0.25 


20 


2 


0.17 


07/22/03 


0.41 


H 


1.66 


0.33 


20 


2 


0.09 


07/21/03 


0.36 


K s 


2.18 


0.35 


20 


2 


0.08 


07/22/03 


0.36 



a DIT corresponds to the exposure time of one image. 
b NDIT corresponds to the number of exposures. 

C PSF corresponds to the average PSF of detected objects in the field. 
d Seeing measured by the ESO Health Check Monitor at each band. 



Table 4. Properties of O-type stars in W51A 



ID 


K s 


H-K s 


A Ks a 


Dzams (kpc) 


D v (kpc) 


Spectral Type 


#44 


11.31 


2.35 


4.08 


1.7 


2.0 


05 


#50 


12.23 


1.88 


3.28 


2.8 


3.6 


06.5 


#57 


10.77 


2.14 


3.72 


1.9 


2.1 


04 


#61 


12.51 


2.65 


4.59 


1.4 


2.0 


07.5 


Average 




2.25 


3.92 


2.0 ± 0.3 b 


2.4 ± 0.4 b 





& A k s was d e termi ned from H — Ks, using [H — Ks)o 
law of lMathisI (ll99ol ) 

b The uncertainty quoted is the error in the mean. 



-0.05 and the extinction 



